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• Premise of the study: Microsatellite loci were developed for the Moroccan endemic endangered species Argania spinosa with 
a combination of a typical library enrichment procedure and a 454 GS FLX Titanium-based high-throughput sequencing 
approach. 

• Methods and Results: A genomic DNA library was enriched and further screened using (GA)i5, (GTA)8, and (TTC)8 biotin- 
labeled probes coupled with chemi-luminescence detection. To increase simple sequence repeat (SSR) loci number, an ultra- 
high-throughput sequencing-based approach was used. Evaluation of all primer pairs was performed with labeled dUTP on an 
ABI 3130x1 sequencer. Eleven polymorphic SSR loci were selected out of 79 SSR regions and extensively characterized on 
150 individuals from eight populations. Total alleles ranged from six to 19 alleles per locus while expected heterozygosity 
ranged from 0.618 to 0.869. 

• Conclusions: The SSRs developed here will be used to further characterize the genetic diversity of A. spinosa across its distri- 
bution range, mainly in the southern part of Morocco and southwestern Algeria. They may also be transferable to other Sapota- 
ceae species. 
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Argania spinosa (L.) Skeels is endemic to Morocco and is 
known worldwide for its oil, which is extremely rich in unsatu- 
rated fatty acids and used in the food and cosmetic industries. 
Argan is the only Sapotaceae species in the region. Argan forest 
is found primarily in the southwestern part of Morocco and in the 
Souss-Massa-Draa region (Msanda et al., 2005); a few popula- 
tions have also been identified in southwestern Algeria (Kaabeche 
et al., 2010). Currently, argan is facing critical regeneration 
problems resulting from overgrazing of young sprouts and trees 
by goats, coupled with inadequate protection; consequently, ar- 
gan tree genetic diversity is severely threatened. Molecular 
analyses of A. spinosa are very recent and still quite limited, 
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mainly due to the lack of species-specific molecular tools. A 
few studies have tried to evaluate its genetic diversity, first with 
the analysis of isozymes (El Mousadik and Petit, 1996a) and 
chloroplast DNA (El Mousadik and Petit, 1996b) and later with 
interspecific random-amplified polymorphic DNA (RAPDs) 
and simple sequence repeats (SSRs) initially developed on 
two different Sapotaceae species (Majourhat et al., 2008). How- 
ever, only six SSRs were polymorphic, and these exhibited 
low diversity according to what would be expected in a peren- 
nial spontaneous (i.e., natural, nondomesticated) tree. Spe- 
cies-specific codominant molecular markers are thus critically 
needed. 

Here, we developed de novo, highly polymorphic and 
A. spinosa-specific SSR markers and used them to evaluate the 
genetic diversity of a collection of geographically diversified 
argan trees. 



METHODS AND RESULTS 

Argania spinosa genomic DNA (gDNA) was extracted from young dried 
leaves and ground using a MM400 mixer mill (Retsch, Dusseldorf, Germany) 
for 2 min at 30 Hz in 1 .5-mL centrifugation tubes with two stainless steel beads 
(Retsch, no. Fr0120). Between 10 and 20 mg of powder was used for gDNA 
extraction using the NucleoSpin Plant II DNA extraction kit according to the 
manufacturer's instructions (Macherey-Nagel, Duren, Germany). DNA con- 
centration was adjusted to 50 ng/jiL. 

Sixteen DNA samples originating from two populations geographically 
separated by more than 400 km were used for library construction (Appendix 1). 
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The protocol of Glenn and Schable (2005) was slightly modified to prepare 
genomic libraries enriched for (GA)i5, (GTA)8, and (TTC)8 repeats. Biotin- 
labeled oligonucleotides were used for library enrichment and Southern analy- 
sis, which were performed with the Biotin Luminescent Detection Kit (no. 
11 811 592 910, Roche Diagnostics GmbH, Mannheim, Germany) and a Bio- 
Rad Fluor-S MAX Multilmager (Bio-Rad Laboratories, Hercules, California, 
USA). Positive colonies containing an insert larger than 400 bp were se- 
quenced by Macrogen (Seoul, Korea) on an ABI 3730x1 genetic analyzer (Ap- 
plied Biosystems, Foster City, California, USA) with BigDye Terminator 
version 3.1 chemistry (Applied Biosystems). SSR motifs were confirmed in 
31 nonredundant sequences. Alternatively, to increase the number of SSR 
loci, gDNA enrichment and direct sequencing using 454 GS FLX Titanium 
(Roche, Basel, Switzerland) were conducted as described earlier (Micheneau 
et al., 2011). 

A total of 79 primer pairs (31 enriched sequences and 48 sequences developed 
through next-generation sequencing) were designed with Primer3 (Rozen and 
Skaletsky, 2000) with the following design parameters: GC content 40-60 (opti- 
mum 50), melting temperature (T^ 55-65°C (optimum 60°C), primer size 18-27 bp 
(optimum 22 bp), GC clamp 1. A typical pigtail (GTTTCTT-) was added to the 
5' end of the forward or reverse primers. Primers were synthesized at Integrated 
DNA Technologies (Leuven, Belgium). They were first tested at 55°C on control 
gDNA. Markers giving no or multiple amplification products were discarded. For 
the remaining markers, each locus was analyzed on a set of geographically distant 
argan trees. For this step, amplification was performed using labeled dUTP nucle- 
otides. A 50-jiL reaction mixture contained 10 jiL of 5x colorless GoTaq reaction 
buffer (Promega Corporation, Madison, Wisconsin, USA), 10 pmol of each primer, 
5 iimol dATP, 5 |imol dCTP, 5 |imol dGTP, 4 |imol dTTP, 0.4 |imol dUTP-DY681 
(Dyomics, Jena, Germany), 1.5 unit GoTaq polymerase (Promega Corporation), 
and 250 ng of gDNA, and was brought to volume with nuclease-free water (Pro- 
mega Corporation). The amplification program was 94°C for 2 min; followed by 35 
cycles at 94°C for 30 s, 55°C for 30 s, and 72°C for 30 s; followed by a final elonga- 
tion step at 72°C for 7 min and then hold at 10°C. Amplification products were 
analyzed on a CEQ 8000 fragment analyzer (Beckman Coulter, Brea, CaUfomia, 
USA), using MapMarker MM-Dl as size standard (Bio Ventures, Murfreesboro, 
Tennessee, USA) and Frag-4 separation method with a separation time of 100 min. 
Amplifications were performed on iCycler (Bio-Rad Laboratories). Out of the first 
79 primer pairs tested, 1 1 fulfilled selection criteria regarding polymorphism, lack 
of background amplification, and ease of scoring. A table summarizing the evalua- 
tion results of all 79 primer pairs is available as Appendix SI. 

Validated SSR markers were resynthesized by Eurofins MWG Operon 
(Ebersberg, Germany) and labeled with 6-FAM, Atto 565, or HEX fluorescent 
dyes. Amplifications were performed in 10-|iL volumes using 0.3 unit of GoTaq 
DNA polymerase (no. M300, Promega Corporation), 2 jiL of 5x GoTaq PCR 
buffer (no. M7921, Promega Corporation), 10-50 ng of gDNA, 0.1 |iL of 
10 mmol/L dNTP mix (no. R0192, Thermo Fisher Scientific, Waltham, Mas- 
sachusetts, USA), and 2 pmol of each primer. ASMS markers (Table 1) were 
amplified individually and further pooled for separation and detection on an 
ABI 3130x1 genetic analyzer (Applied Biosystems), while ASMS2012 markers 
were multiplexed for amplification in two triplex reactions. All amplifications 
were performed according to the melting temperatures reported in Table 1 . 

One hundred fifty argan trees originating from seven populations located in 
southern Moroccan argan forests and one population from the Oued Grou relict 
area, in northwestern Morocco, were used for in-depth testing of the newly de- 
veloped SSRs (Appendix 1). Allele number (A), allelic richness (A^), and ex- 
pected (//e) and observed {H^ heterozygosities were calculated using FSTAT 
software version 2.9.3.2 (Goudet, 1995). Results revealed high allele numbers 
(Table 1). A total of 136 alleles were scored. A ranged from six (ASMS2012-10) 
to 18 (ASMSOl). A, varied from 3.992 to 9.085. and ranged from 0.618 
to 0.869 and 0.532 to 0.887, respectively. No significant difference was observed 
between and for nine loci, indicating that populations are approaching 
Hardy- Weinberg equilibrium. A difference has been noticed for loci ASMS 19 



and ASMS 3 1 , where was higher than H^, indicating a slight lack of heterozy- 
gotes. This might be due to natural selection of loci near the two SSRs or to the 
presence of null afleles. 

Argan tree is considered as a spontaneous natural species with no genetic 
change due to human selection. The observed levels of genetic diversity are thus 
explained by the combined effects of mutation, genetic drift, gene flow, and 
natural selection. 



CONCLUSIONS 

In our study, 150 argan trees were genotyped with 11 SSR 
markers. Samples originated from eight geographically sepa- 
rated regions with no possibility of natural gene exchange. The 
1 1 microsatellite markers used in this study were very informa- 
tive and showed a high ability to detect polymorphism within 
analyzed genotypes. Outcomes of this research provide the first 
codominant markers specific to A. spinosa and will help in elab- 
orating strategies to better protect this threatened species. They 
will contribute more efficiently to evaluate the genetic diversity 
and to understand gene flow among and between argan forest 
areas. Moreover, developed markers may constitute valuable 
tools for evaluation of genetic diversity in other Sapotaceae 
species. 
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Appendix 1. Location information for Moroccan populations of Argania spinosa used in this study/ 



Population 


Site name 


Region 




Latitude (°N) 


Longitude (°W) 


Altitude (m) 


POPOl* 


AGD 


Agadir 


23 


30°24'56.3" 


009°36'17.3" 


32 


POP02* 


OGl 


Rabat cued grou 


29 


33°55'48.2" 


006°45'00.0" 


429 


POP03 


ESSII 


Essaouira 


20 


31°55'35.3" 


009°24'28.9" 


108 


POP04 


AGIOSI 


Agadir 


20 


30°3ri4.1" 


009°4r42.3" 


34 


POP05 


AG10S5 


Agadir 


21 


29°0r54.0" 


009°58'07.6" 


412 


POP06 


B2-02 


Agadir 


17 


30nr52.5" 


009°18'53.5" 


178 


POP07 


B3-09 


Agadir 


12 


30°04'28.0" 


009°09'07.4" 


563 


POP08 


EI 


Essaouira 


8 


31°24'58.8" 


009°42'58.5" 


97 



Note: N= sample size. 

^No herbarium vouchers were collected due to the sensitive status of the species. 
* Population used for DNA library generation. 
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